Quasi-exactly solvable quartic Bose Hamiltonians 

S. N. Dolya 1 and 0. B. Zaslavskii 2 
l B. Verkin Institute for Low Temperature Physics and 
Engineering, 
47 Lenin Prospekt, Kharkov 61103, Ukraine 
E-mail: dolya@ilt.kharkov.ua 
2 Department of Mechanics and Mathematics, Kharkov V.N. Karazin's 

National 
University, Svoboda 
Sq.4, Kharkov 61077, Ukraine 

E-mail: aptm@kharkov.ua 
We consider Hamiltonians, which are even polynomials of the forth or- 
der with the respect to Bose operators. We find subspaces, preserved by 
the action of Hamiltonian These subspaces, being finite-dimensional, include, 
nonetheless, states with an infinite number of quasi-particles, corresponding 
to the original Bose operators. The basis functions look rather simple in 
the coherent state representation and are expressed in terms of the degen- 
erate hypergeometric function with respect to the complex variable labeling 
the representation. In some particular degenerate cases they turn (up to the 
power factor) into the trigonometric or hyperbolic functions, Bessel functions 
or combinations of the exponent and Hermite polynomials. We find explicitly 
the relationship between coefficients at different powers of Bose operators that 

ensure quasi-exact solvability of Hamiltonian. 
PACS numbers: 03.65.Fd, 03.65.Ge 



I. INTRODUCTION 



The conception of quasi-exactly solvable (QES) systems, discovered in 1980s M - [15 



received in recent years much attention both from viewpoint of physical applications and 
their inner mathematical beauty. In turned out that in quantum mechanics there exists a 
peculiar class of objects that occupy an intermediate place between exactly solvable and 
non-solvable models in the sense that in an infinite Hilbert space a finite part of a spectrum 
is singled out within which eigenvectors and eigenvalues can be found from an algebraic 
equation of a finite degree - in other words, a partial algebraization of the spectrum occurs. 
For one-dimensional QES models corresponding QES Hamiltonians possess hidden group 
structure based on sl(2, R) algebra. Thus, they have direct physical meaning, being related 
to quantum spin systems |13| . 

Meanwhile, the notion of QES systems is not constrained by potential models and can 
have nontrivial meaning for any kind of infinite-dimensional systems. In the first place, it 
concerns Bose Hamiltonians whose physical importance is beyond doubt. Here one should 
distinguish two cases. First, it turns out that some systems of two interacting particles 
or quasi-particles with Bose operators of creation and annihilation a, a + and b, 6 + can be 
mapped on the problem for a particle moving in a certain type of one-dimensional potentials 
and, remarkably, these potentials belong just to the QES type [0] - |TB|. In particular, such a 
type of Hamiltonians is widely spread in quantum optics and physics of magnetism |13| . The 
aforementioned mapping works only for a special class of Bose Hamiltonians which possess 
an integral of motion. Then the procedure is performed in three steps: (i) all Hilbert space 
splits in a natural way to different pieces with respect to the values of an integral of motion, 
(ii) in each piece the Schrodinger equation takes a finite-difference form, (iii) it is transformed 
into the differential equation by means of introducing a generating function. In so doing, the 
integral of motion under discussion represents a linear combination of numbers of particles 
a + a and b + b. 

The second kind of Bose systems looks much more usual - it is simply some polynomial 



2 



with respect to Bose operators of creation and annihilation of one particle. The fact that 
only one pair a, a + enters Hamiltonian, deprives us, by contrast with the first case, of the 
possibility to construct a simple integral of motion - in this sense the eigenvalue problem 
becomes more complicated. In general, the solutions of the Shrodinger equation contain in- 
finite numbers of quasi-particles and only approximate or numerical methods can be applied 



to such systems. However, as was shown recently |L7|], if the coefficients at different powers 
of a, a + are selected in a proper way, in some cases a finite- dimensional closed subspace is 
singled out and algebraization of the spectrum occurs similar to what happens in "usual" 
QES potential models or differential equations. In so doing, the eigenvectors belonging to 
the subspace under discussion, can be expressed as a finite linear combination of eigenvectors 



of an harmonic oscillator and, thus, contain a finite number of quasi-particles [17 



In the present article we extend the approach of |17] and consider much more general 
classes of Hamiltonians. Their distinctive feature consists in that the relevant basis functions 
that compose a finite-dimensional subspace, look very much unlike the wave functions of a 
harmonic oscillator. As a result, we obtain QES models with an infinite numbers of quasi- 
particles in this /imte-dimensional subspace. Bearing in mind physical application, we make 
emphasis on Hermitian Hamiltonians, although our approach is applicable to more general 
QES Bose operators without demand of Hermiticity. 

II. BOSE HAMILTONIANS AS DIFFERENTIAL OPERATORS AND 
STRUCTURE OF INVARIANT SUBSPACES 

Consider the operator which is the even polynomial of the forth degree with respect to 
Bose operators of creation a + and annihilation a. It can be written in the form 

H = a ++ K 2 + + a__K 2 _ + a 00 K 2 + a _K K^ + a +0 K + K + a K + a„K„ + a + K + , (1) 

where 

1 / , 1\ a 2 a+ 2 



[K , K ± ] = ±K ± , [K+, K„] = -2K . (3) 

The Casimir operator C = K I - \(K + K^ + K_K+) = 
We will use the coherent state representation in which 

a ~* Yz a+ ~* z (4) 

After substitution into (|l|) Hamiltonian H(a + ,a) becomes a differential operator H(z,-§z). 



In the previous article |17| we discussed Bose systems that possess the invariant subspace of 
the form F = span{z n } or span{z 2n }. The first natural step towards generalization consists 
in considering subspaces (with N fixed ) 

F = span{u n }, u n = z 2n u, n = 0, 1, 2...N, (5) 

for which the following procedure should be realized, (i) The action of operators of Ki on 
the functions u n should lead to the linear combinations of functions from the same set {u n }, 
(ii) by the selection of appropriate coefficients in ([I]), we achieve the subspace F to be closed 
under the action of Hamiltonian H. We would like to stress that the condition (i) does 
not forbid u n with n > N to appear in terms like KiU n but the condition (ii) rules out 
such functions from Hu n (recall that we consider Hamiltonians which are quadratic-linear 
combinations of K^). 

It is seen from (Q), (^j) that the operators Ki contain z and Therefore, it is convenient 
to assume that differentiation of u(z) gives rise to u up to the factor that contains powers 
of z. The corresponding choice is not unique. In the present article we restrict ourselves to 
one of the simplest possibilities that leads to nontrivial solutions. To this end, we choose u 
that obeys the differential equation 

u' = A(z)u, A(z) = (- + 2pz). (6) 

z 

We will show below that the choice (|) relates KiU n to u n , u n ±i that, in turn, allows us to 
formulate the conditions of cut off for Hamiltonian in the form of algebraic equations which 
its coefficients obey. It follows from (g) that u = z 13 exp(pz 2 ) . To ensure asymptotic analytic 



behavior near z = 0, we demand that (3 = 0, 1, 2... Now let us take into account some basic 
properties of coherent states (see, e.g. Ch. 7 of Ref. Our functions tt n (z) must belong 

to the Bargmann-Fock space. It means that they should obey the conditions of integrability 
and analyticity. The condition of integrability for any two functions /, g from our space 

dzdz*f*ge~ zz * < oo, (7) 

entails, for our choice of u, \p\ < 1/2. 

Taking into account eq. (^), it is straightforward to show that 

K + u n = C + u n+ i, (8) 
K-U n = A_ (n)u n + B_ (n)w n _! + C_w n+1 , (9) 
K u n = A (n)u n + C u n+ x, (10) 



where C + = §, C_ = 2p 2 , C = p, A_(n) = (2/5 + 4n + l)p, A (n) = 2/3+ f +1 , S_i 



(/3+2n)(/3+2n-l) 
2 

Using eqs.(|S|) - (|T0"D, one can present the action of the operator (TJ) in the form 

Hu n = D 2 u n+2 + Di(n)u n+ i + D (n)u n + Di(n)w n _i + D 2 (n)u n - 2 , (11) 

where 

^2 = + a -C_C + a 00 C 2 + a__C_C_ + ^±± (12) 

Di (n) = a__ [A_(n)C_ + C_A_(n + 1)] + a o [4)(7i)C + C A (n + 1)] + ^A (n) + (13) 
a _[A_(n)C + C_A (n + 1)] + a C + a_C_ + 

Do = a 00 A 2 (n) + a -[A^(n)A (n) + B_(n)C ] + a A (n) + a^A^(n), (14) 
D l (n) = a__ \A_(n)B_ (n) + B_ (n) A_(n - 1)] + a „B_(n)A (n - 1) + a_5_(n), (15) 

D 2 = a__B_ {n) B_ (n - 1) . (16) 

For the operator ([!]) to be quasi-exactly solvable with the invariant subspace (^), it is 
necessary that the following conditions of cut off be satisfied: 



D 2 = 0, 
D 1 {N) = 0, 
#i(0) =0, 
Aj(0) = o, 
D a (l) =0. 



(17) 
(18) 
(19) 
(20) 
(21) 



In general, this system is rather cumbersome. However, it is simplified greatly if we consider 
Hermitian Hamiltonians with a = = a ++ . Then D 2 = and we get three equations 



p 2 a 00 + 2p 3 a - + ^pa - = 



r «o- 



■(2/3-3)+a_]/3(/3- 1) = 



(22) 
(23) 



1+V r 2/? + 4iV + l 2 7 + 6/3+12JV. a 00 

2 — a- + pa + [— + p 2 ^ ]a _ + -yp(2/3 + AN + 3) = 0. (24) 

It is assumed that all coefficients are real. The analysis leads to the 

following table of possible solutions: 
a ++ = 





P 


P 




&00 


a 


1 








4 a 0- 


a.v. 


a.v. 


2 





1 


47V+3 _ 
4 a 0- 


a.v. 


a.v. 


3 


a.v. 


a.v. 


3-2/3 
4 a 0- 


f(p)ao- 


[Wjr+i 2(iV + i) p]ao _ 



4 


a.v. 





a.v. 


f(p)ao- 


/(p)o_ + /i(p, iV)a _ + / 2 (p, iV)a ++ 


5 


a.v. 


1 


a.v. 


f(p)ao- 


/(p)a_ + /i(p, iV + ±)a _ + / 2 (p, AT + i)a ++ 



Here "a.v." denotes "arbitrary value" with the reservation that /3 is a positive integer or 
zero and |p| < 1/2, as is explained above. By definition, 

4iV + 3 



/(p) 



l+4p 2 4iV + 5 4iV + l 

fi[P,N) = — — p, f 2 (p,N) 



■(16p 4 -l). (25) 



2p ' J vr ' 7 8p 2 »p- 

One can check that the solution 1, when u — 1, corresponds to even states for the 
example considered in eq. (7) of [|17|] provided in that equation the coefficient A 2 = 0. In a 



similar way, the case 2 (u = z) corresponds to odd states from the same example. However, 
the cases (3)-(5) represent new solutions that were not contained in JET]. 



III. DOUBLED INVARIANT SUBSPACES 

Much more rich family of new classes of Bose QES quartic Hamiltonians can be obtained 
if we generalize the structure of the invariant subspace introducing, in addition to «„, a 
second subset of independent functions. Consider the set of functions 

u n = z 2n u, v n = z 2n+1 v, n = 0, 1, 2... (26) 
We are interested in such function which form a set, defined for a fixed N = 0,1, 2..., 



F = span {u n , v n } = span {z 2n ■ u (z) , z 2n+1 ■ v (z)} , (27) 
n = 0,l,2,...,N ; 



invariant with respect to the action of the operator H. The dimension of F is equal to 
2(N + 1). It may happen that, for some values of parameters, the functions v n may be 
proportional or even exactly equal to u n . Then our subspace reduces to the N + 1 one (|5[) 
considered in a previous section. To gain qualitatively new QES models, in what follows we 
will consider the functions u n and v n as, generally speaking, independent. 

Let u and v obey the system of differential equations that generalizes the relation (^|): 

u = Au + Bv, (28) 
v' = Cu + Dv. (29) 

Here prime denotes differentiation with respect to z, A, B, C, D are functions of z. 
It follows from (BSD, (E§D that 



(30) 



✓ -✓[*4W + M] = 0. (31) 



Here S = A + D = SpL, where L 



\ CD ; 



, W(fi,f 2 ) = f[f 2 - /1/2 is a Wronskian, 



A = AD — BC is a determinant of L. In what follows we assume for simplicity that quantities 
in denominators in fl30l), (|3ll) B(z) = a = const and C(^) = 7 = const. Then we have 

«" - u'S + u(A - A') = 0, (32) 
v" - v'S + v(A - D') = 0. (33) 

We assume also, by analogy with (|B|), that A and D contain only terms of the order z and 
z~ x in the Loran series: A = 2pz + f3z~ l , D = 2tz + 5z~ x . Then we have 

d (3 

—u (z) = av (z) H — w (2) + 2pz« (2) , (34) 
dz z 

d S 
—v (z) = 7M (z) H — v (z) + 2tzv (z) . 

CtZ z 

One obtains from ((53), fl53]) 
it" - u'[^±^ + 2(p + t)z) + m ^ 1 ^"^ + 4prz 2 + 2f3r + 2p5 - 2p - aj] = 0, (35) 

Z Z 

v " _ v '[*±£ + 2(p + t)z\ + uf ^t + 4prz 2 + 2r/3 + 2p5 - 2r - try] = 0. (36) 

JZ z 2 

Our functions u n (^) and v n (z) must belong to the Bargmann-Fock space that entails, simi- 
larly to what is obtained in the previous section, the conditions |p| < 1/2, \t\ < 1/2. 

To elucidate what constraints are imposed by the demand of analyticity, consider sep- 
arately several different cases. If 7 = a = 0, eqs. ( |34"|) can be integrated and one easily 
finds that u = z 13 exp(pz 2 ), v = z s exp(rz 2 ), whence it is obvious that (3 = 0,1,2... and 
5 = —1, 0, 1, 2... If a = but 7 7^ 0, one can make the substitution v — z s exp(rz 2 )w. Then 

w' = -fZ (3 ~ 5 exp[(p-r)z 2 ] (37) 

It is clear that (3 = 0, 1, 2..., whereas 5 is arbitrary except 5 — (3 + 1 since the latter would 
have led to the logarithmic terms in v(z). The similar situation occurs when 7 = but 
Then 5 = —1, 0, 1, 2... and forbidden values of (3 are (3 = 5 + 1. 
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Let now 07 7^ 0. First, consider the case p 7^ r. Then by substitutions 

u(z) = y [z 2 (r — p)) • exp (pz 2 ) ■ z$ 



(38) 



v(z) = y (z 2 (p - t)) ■ exp (rz 2 ) ■ z s (39) 

eqs. (|35D, (p6|) are reduced to the form, typical of a degenerate hypergeometric function 

d 2 d 
X dx 2 ~ V ^ + ( T] ~ x ) -^v ( x ) - & ( x ) = °> ( 40 ) 

where rj — ^(/3 — 8+1),^ — 4( yZ p ) ■ The function ?/ satisfies the equation of the same form 
(^0|) but with parameters f]= 1 — 77, £ = — £. 

To determine the admissible range of parameters (3, 5 one can appeal directly to the 
well-known properties of this function and take into account that the general solution of eq. 
(flOl) has the form y = Ay x + By 2 , where y x = $(£, r/; x) and y 2 = x 1_r? $(£ — rj + 1, 2 — rj; x) 
and the standard notation for the degenerate hypergeometric function is used (see Ch. 6 
of Ref. [0). First consider the case when t] is non- integer. Then $ — > 1 when x — > 
and from fl38|) we obtain the function u can have two possible asymptotic forms: u\ ~ z 13 
and U2 ~ z +1 . The function f (2) behaves, correspondingly, like t>i ~ 2; /3+1 and v 2 ~ Z* 5 . 
Therefore, it turns out that there are two cases: 

/3 = 0, 1, 2..., 5 is arbitrary (41) 



or 



5 = -1, 0, 2...,P is arbitrary. (42) 

If r] is integer, there exists only one independent solution of eq. (|40|), regular at x — > 0. 
The corresponding solution is known to be $*(£,?7;x) = rfr?f^ ' ^ n ^ ne li m ^ x — * 
$* ~ x 1 ^ 71 that does not affect the conclusion about admissible range of j3 and S. 

Let now p = r, 0:7 7^ 0. By substitution 

w(,2) = y (2^/aPyz) ■ exp (z (pz — \forf)) ■ z@ 



we obtain that the function y(x) obeys the equation 

d 2 , x , n , N d , s (/? — 8) . . 

x i~^y ( x ) + (p-°- x ) i~y ( x ) — o — v w = 

ax z ax 2 

that has the same form as (HO) and admissible /3 and 5 satisfy one of criteria ([II]), (f^)- 

Differential equations for our functions can be also written in the symmetric form. Let 
us make the substitution 



Then 



where 



T §+s r (p + r) 2l 
u = 2 exp[— ]. 



+ (e - K//)* = o, 



^e// = ^ + (p-r)V, 
(5-/5)(5-/5 + 2) 



(43) 
(44) 

(45) 
(46) 



e = (p — t)(5 — (3 — 1) — cry. 



(47) 



In a similar way, 



~ g+l Ap + r) 2 



(48) 



where ^ obeys the equation ( [44]) with the same structure of V e ff but with another k 



(P-8){(3-6+2) 



(r — p){(3 — 5 — 1) — cry. It is seen that e can be obtained from e by 



interchange between p and r, j3 and <5. It is also seen that e — e = 2{p — r). 

Thus, formally, we obtain the harmonic oscillator with a barrier z~ 2 (Kratzer Hamilto- 
nian). Let us remind, however, that the variable z in our context is complex. 

It follows from ([34]) and ([43]), (f48|) that the functions ^ obey the system of equations 

(49) 



d , s 8 — 3 



d , , 3 — 8 
5 7 + (p-r), + — 



(50) 
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IV. CONDITIONS OF QUASI-EXACT SOLVABILITY 



The action of operators K { in the subspace ( P7|) has the same structure (|J) - fllCf ) but 
now the corresponding quantities A iy B i: Ci becomes 2x2 matrices: 



K+fn — C+Zn+lj 



(51) 



K^f n = A^{n)f n + B_(n)/ n _i + C_/„ +1 , 



(52) 



ffo/n = 4W/n + Co/n+l, 



where /„ 



c 4 



V 01 / 



2 



2p 2 
\j(r + p) 2r 2 y 



C = x 



A_(n) 



/ 



2p 

7 2r 

\ 



4/5p + 8np + 2p + «7 2a(r + p) 
^ 7 (/3 + 5 + 2 + 4n) 45r + 8nr + 6r + cry y 



/ 



2(3 + An + 1 2a 

25 + 4n + 3 



/ 
\ 



/ 



(53) 



(54) 



(55) 



(56) 



(57) 



(P + 2n)((3 + 2n - 1) a(f3 + 5 + An) 
(5 + 2n + l)(<S + 2n) 

In a similar way, the action of Hamiltonian in the invariant subspace can be represented 
in the form 



Hf n = D 2 f n+2 + D l {n)f n+l + D (n)f n + Di(n)/ n _! + D 2 (n)f n „ 2 , 



(58) 
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where matrices Di and Di have the form (|l"2"D - ( |I~6"D with Aj, taken from eqs. fl54|) - 

(|57|) (the order of operators is taken into account properly in this form of writing). 

For the operator (|l|) to be quasi-exactly solvable with the invariant subspace ([27]), it is 
necessary that the matrix version of the conditions of cut off ([[7]) - fl2~T|) be satisfied. The 
corresponding system of equations is too cumbersome to be listed here. It can be simplified 

greatly if we assume the condition a = 0, in which case after simple calculations we get 

D 2 = 0, 



( 



\ 



p 2 a 00 + 2p 3 a - + 7;pa + o + ^p 
y lY T 2 a 00 + 2r 3 a _ + §ra +0 + pp J 

where Y = [~(p + r)a o + (p 2 + t 2 + pr)a - + \a +0 
( 

2D X (0) = 



(59) 



\O-0- 



(2/3 - 3) + a_]/?(/3 - 1) a[^(2(3 2 + 25 2 + 2(35 - 3(3 - 5) + a_(/3 + 5)] 







<5(<5+l) 



[(2tf-l)5^ + a 



(60) 



£>i(iV) 



\ <^3 ^4 y 



(61) 



a+ 



+ 2p 2 a_ + pa + ^(2/3 + 4iV + 1) + ^-[a^r + 2p) + p 2 (7 + 6/3 + 12JV)] (62) 



«o 



+^[p(4/3 + 8iV + 6) + a7], 



(63) 



do = aY, 



(64) 



4 = Z = [2(t + p)a_ + a + ^£ + ^(/3 + 5 + 4iV + 4)], (65) 



£ = 12iV(p + r) + a 7 + r(45 + 2(3 + 11) + p(4/3 + 25 + 9), 



(66) 



4 = -it + 2r 2 a_ + ra + -p(25 + 4iV + 3) + -^[arf(p + 2r) + r 2 (12iV + 13 + 65)] (67) 
2 8 2 



+ ^[r(45 + 8iV+10) + a7] 



(68) 
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One can observe that 

d x - d 4 = - t)Z + ((3 - 5 - l)Y (69) 

The system of equations ( |iT|) - (|21| ) with ([^) - ( pTf ) taken into account looks rather cum- 
bersome but the relation ( |B"9"D simplifies analysis significantly. 

However, for a generic case a ^ algebraic calculation are so bulky that we had to 

resort to using a computer. 

In what follows we restrict ourselves by the Hermitian case only which is the most 
interesting for physical applications. This implies that 
where all coefficients are real. 

The full set of non-trivial Hermitian solutions of the system ( |T7D - fl21|) and their classi- 
fication are given in the next Section. 

V. HERMITIAN SOLUTIONS OF ALGEBRAIC EQUATIONS 

For the type of invariant subspaces (^) under consideration, we suggest below the clas- 
sification of all QES Hermitian Hamiltonians, quadratic-linear with respect to operators Ki 
(i.e., those which represent even polynomial of the fourth order in terms of a, a + ). 

In the tables below we list only qualitatively different cases in the following sense. If 
some solutions can be obtained by the limiting transition (a ++ —>■ 0, 7 — > 0, etc.) from 
a more general case, we do not repeat them. As before, we use abbreviation "a. v." for 
" arbitrary value" . 

A. a ++ = 

7^0 



13 





5 


P 


a 


P 


7" 


a 


Onn 


fin 


1 


a.v. 





a.v. 


o 


cey 


1-25 On 


/ (t) «o- 




2(2/V+l-5) 


4 "0— 


I'yl"/ T 4r J^O- 


2 


a.v. 





a.v. 









/ (p) a o~ 


[/, (g) + <5+4JV+4i 


4(AT+1) 


4 t *0— 


Lpy v u / ~ 4p j u o- 


3 


a.v. 


1 


a.v. 


o 


07 




f (t) a _ 


[ra C<5 1) + 5+4Ar+5 lan 


2(2JV+2-<5) 


4 "0- 


L' y \ u L ) ~ 4 T j^o- 


4 


a.v. 


1 


a.v. 







l-2<5 


/ (p) a o- 




4(iV+l) 


4 a 0- 


[py (* - 1) + 


5 





a.v. 


a.v. 





ay 


3-2/3 


f (r) a _ 


[Tff(/? 2) + ^+ 4 ]a _ 


4(AT+1) 


4 a 0- 




6 





a.v. 


a.v. 







3-2/3 


/ (p) ao- 




2(27V+3-/3) 


4 a - 


[p^(/3-2) + ^± 4 ]a _ 


7 








a.v. 





a.v. 


a7-r(4Af+l) 


/ (r) a _ 


(-r-a 7 + ^)«o- 


8 








a.v. 


a. v. 





ay-p(4,N+3) „ 

4p °0- 


/ (p) a o- 


(p - a 7 + ^±i)a _ 


9* 








a.v. 








— — a o 


a. v. 


-2(iV+l)aoo 


10 





1 








a. v. 


4/V+3 _ 

4 a o- 


/ (t) a - 


4/V+5 _ 
4r °0- 


11 





1 





a.v. 





4AT+3„ 
4 a 0- 


/ (p) a o- 


4/V+5 
4p fl 0- 


12 


-1 


1 


a.v. 





a.v. 


a7-r(4Af+3) 

4r fl 0- 


f (t) <»o- 


(r 07+^)00- 


13 


-1 


1 


a.v. 


a.v. 





a 7 -p(47V+l) „ 

4p °0- 


/ (p) oo- 


(-p - « 7 + ^)ao- 


14 


-1 


1 


a.v. 








-^ann 


a. v. 


— 2(N + lWn 


15 


-1 








o 


a.v. 


4 "(J— 


/ (t) <»o- 


4r "°- 


16 


-1 








a.v. 





4AT+1 
4 «0- 


/ (p) oo- 




17 


-1 


a.v. 


a.v. 





ay 


3-2/3 


/ (t) «o- 


[ry 09 - 1) + ^P]«o- 


4(AT+1) 


4 oo- 




18 


-1 


a.v. 


a.v. 







3-2/3 


/ (p) ao- 




2(2AT+2-/3) 


4 a 0- 


[p^(/3-l) + ^K- 



where 



g ( x ) = x - 3 - AN (70) 

and we used definition of / according to (§5|). 

*Note: in the case 9 the coefficients a +0 = a _ = 0. 
7 = 
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5 


P 


a 


P 


T 


a_ 




a 


19 


a.v. 


5 + 1 





a.v. 


P 


1-25 
4 a O- 


/ (p) o - 




20 





1 











(47V+3) 
4 °0- 


a.v. 


a.v. 


21 


-1 














(1+4JV) 
4 a O- 









B. a ++ ^ 

In all admissible cases 22 - 25 p and r take arbitrary values, 

o - = -^-(r + p) (4rp + l)a ++ ,a 00 = ((4 r P + l) 2 + 4 ( p 2 + r 2 ))a++. 



The rest of relevant quantities is 





5 


& 


a 


a 


a_ 


22 








a.v. 


MP, r ) a + + 


/_(p,r)a ++ 


23 





1 





g (p,T, N)a + + 


#_(p, r,iV)a + + 


24 


-1 


1 


a.v. 


fo(p, T ) a + + 


/_(r, p)a ++ 


25 


-1 








9o(p,t,N- \)a + + 


#_(p, r,iV- ±)a+ + 



Here / (p, r) = ^[(16r 2 p 2 - 1) (N + 1) - p 2 + r 2 + a 7 (r + p)], 

Sb(p, r, iV) = &*p. (16r 2 p 2 - 1), (p, r, iV) = [rp(lQN + 28) - 4iV - 3] 

/_ (p, r) = --L. [r 2 p(16iV + 28) + r p 2 (20 + 16iV) + a 7 (4pr + 1) - p(4iV + 3) —r(l + 4N)] 

7 = 

Now p = r. 





5 





a 




a _ 


a 


26 








a.v. 


/oo (r)a+ + 


2/(r)a + + 


-/ 2 (r)a ++ + /(r)a_ 


27 





1 





/(r)a _ + /(2r 2 )a ++ 


a.v. 


/(r)o_ - # (t, r, iV)a ++ + /i(r, JV + |)a _ 


28 


-1 


1 


a.v. 


foo(r)a + + 


2/(r)a+ + 


-/ 2 (r)a ++ + /(r)a_ 


29 


-1 








/(r)ao_ + /(2r 2 )a ++ 


a.v. 


/(T)a_ - g (r, T,N— |)a ++ + /i(r, JV)a _ 



/oo(r) = -^(l + 16r 4 + 16r 2 ), 



the function fi is defined according to (|25|) . 
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VI. EXPLICIT EXAMPLES OF INVARIANT SUBSPACES 



In this section we list shortly the explicit form of solutions for u(z) and v(z). 

1) The case 9: P = 5 = p = r = 0, 

a) 7 = —a = —uj 7^ 

Then it follows directly from that the functions u, v can be chosen as u = cosujz, v = 
sina;2;. After simple calculations one finds that, apart from the Hermitian QES Hamiltonian, 
there exists also the non-Hermitian QES operator: 

H = a 00 K 2 + a -K K^ + \^-lu 2 - 2a 00 (N + 1)]K + a_K_ + ^Vtf+ (71) 

b) In a similar way, we obtain for 7 = a = uj that u = cosh (uz), v = sinhu;;?, and the 
operator H is obtained by replacement uj 2 — > —uj 2 in the expression ([Tj]). 

2) p = T = 0, 8 = -1 - n, (3 = n (n = 0,1...), a = -1, 7 = 1 

Then we have the following solutions of ([34]) : u = J n (z), v(z) = J n+ x(z) (Bessel func- 
tions). 

(3 + 2n) a 00 a _ a oo, AAT . q n / vo n 
a_ = a -, a +0 = = a++, a+ = — , a = — ^~( 4iV + 3 )- ( 72 ) 

The effective Hamiltonian is non-Hermitian. 

3) Consider the case 8 of Hermitian Hamiltonians (the case 7 can be considered in 
a similar manner): (3 = 5 = = r, then k = = k and formally we have in the z- 
representation the wave function that looks like that of a pure harmonic oscillator would 
look in the coordinate representation. Eqs. ( liSf) , ( j50|) take the form 

fy ~ 

= -j=V, (73) 



2p 
27 



6 + # = -^=* , (74) 



where b — -^[^ - pz], 6 + = -^[J; + P^]- It is obvious that [b, 6+] = 1. 

The frequency is equal to uj = 2p, e = —^ — 0:7, e — e = 2p. Let also e = e n = u(n+l/2), 
= —uj(n + 1). Then we have the eigenvalue and e = e n+ i. Thus, our subspace is 
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span{^ n z 2n exp[(|) z 2 ], ^ n+ iZ 2n+1 exp[(|) z 2 ]}, where \l/ n is the wave function of the n-th 
level of the harmonic oscillator, \l/ n = exp(— ^pz 2 )H n (z^/p), H n is the Hermite polynomial. 
Here (5 = 0, 1, ... 

We would like to stress that our system represents an anharmonic (not harmonic!) Bose 
oscillator. The functions, which have the same form as those of an harmonic oscillator, 
appear in this context in the coherent state representation (not in the coordinate one, as 
would be the case for the usual harmonic oscillator) and represent auxiliary quantities. 



VII. GENERALIZATIONS 

In these section we describe shortly, on the basis of the suggested approach, some pos- 
sible ways of generation of new invariant subspaces, suitable for constructing QES Bose 
Hamiltonians. As the method of constructing is the same as was used above, we only dwell 
upon the structure of subspaces. 

1) Let us introduce quantities 



fi = ^ 2 v, f n = z 2 ^\v, f n = , 2 v 



(75) 



and consider the subspace Fn = span{f^, f 2 , f%}, n = 0,1, ...,N ; N = 1,2,.... Let /„ 
{ f i\ 

J n 



■J n 

c3 



. Then (^Tj)-(^) take place, where, however, now the corresponding matrices have 



dimension 3x3: 



A (n) 



I is a unit matrix, 



U+hot ^ 



q zg-t-z?7-t-a q 



2g+2?j+3 

7 V+jJ 



Co 



A_(n) 



2p (4f + 1) + a 7 2a (u + 2p) 



^ 2p X 



7 , , a 
2 U 2 



2t 



a 



\ 



(2 + 3£ + rj) to (2f + 2r] + 3) + 2a 7 § (2 + f + 3r?) 



2 7 (u + 2t) 
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2t (1 + 477) + cry , 



(76) 



(77) 



Bin] 



( 



£(2£-l) a(3Z + V ) ^ 

7(3?? + 77 ( 277 — 1) J 



(78) 



^8p 2 





7(cu + 2p) 2cj 2 a(2r + cu) 
8r 2 



\ 



£ = /3 + n, r] = n + 5, u = p + r. 

2) Further generalization of consists in considering 



(79) 



fl = z 2n uu, f 2 n = z 2n+1 uv, f* = z^uv, z^uv, f n = z ln vv 



2n + 1 o,^ ^2n+l~ 



~2n ~ 



(80) 



Here u, v refer to the functions obeying eqs. ( p4|) with parameters a, (3, 7, 5, p, f. Then we 
can construct span f 2 , ft fn) ■ 

3) Consider span{z 2n+ ^ kS> u M ~ k v k }, £(k) — 1, if k is odd and £ = 0, if k is even. Here 



fc = 0, 1,2...M, n = 0,l..JV. 



The functions -{V} that appear in the subspaces l)-3) obey the system of equations of 
the type ^ = C) (z) ui that contains, el. S 9 particular case, eqs. 



VIII. SUMMARY 

Let us summarize the basic features of our approach to constructing invariant subspaces 
for Bose Hamiltonians. We consider systems whose Hamiltonian can be expressed in terms 
of generators of Ki (]2|). Further, we (1) use the coherent state representation in which all 
Bose operators become differential ones, (2) split the space of states of an harmonic oscillator 
to even and odd states, (3) deform each of two pieces by introducing, as a factor, additional 
unknown function which is different for each piece, (4) demand that both these functions 
obey a coupled system of differential equations that the action of Ki convert each of basis 
vectors into a linear combination of vectors of the same type, (5) select coefficients, with 
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which Ki enter Hamiltonian H, to ensure the cut off in the space of basis functions. In a 
sense, we introduce a kind of an additional degree of freedom - effective "spin" s. Then 
kinds of subspaces considered in our paper can be assigned the values s = (Sec. II), 1/2 
(Sec. Ill), 1 (Sec. VII), etc. However, we want to stress that this "spin", in contrast to 
matrix QES models |2(J - |23|] , does not appear in Hamiltonian and serves to describe the 
structure of solutions only. 

It is worth stressing that it is just H itself is quasi-exactly solvable, whereas generators 
Ki themselves, with the help of which H is built up, are, in general, not. This is in a sharp 
contrast with "usual" QES Hamiltonians in quantum mechanics. Let us remind ourselves 
that in the latter case H can be expressed in terms of operators of an effective spin Si that 
realize sl(2, R) algebra, each of them possessing finite-dimensional subspace. 

The essential ingredient of our approach is using the coherent state representation. For- 
mally, formulas (|J) look like those in the coordinate-momentum representation. However, 
for our Hamiltonian one should check carefully the normalizability of solutions that im- 
plies integration over whole complex plane in a scalar product and impose constraints on 
admissible values of parameters. Another constraints stem from the demand of analyticity. 

We want to point out that non-Hermitian QES operators can also be of interest. They 
may be used, as auxiliary quantities, in physical applications for finding spectra of Hermite 
Hamiltonians. For instance, they may appear in mappings like HL = LH', where H is 
Hermitian. Knowing the spectrum of H' or its part due to its quasi-exact solvability, one 
can restore a part of spectrum of a physical Hamiltonian H. Apart from this, the approach 
considered in the present article opens a way to the search and classification of linear dif- 
ferential operators with different invariant subspaces. This would enable one to generalize 



or extend the results obtained for such subspaces with a basis of monomials ||24|| , p5| . In 
particular, for the case ([7^) we obtained solutions in the form of combinations of Bessel 
functions and monomials. 

In this paper we restricted ourselves to one-particle systems but the suggested approach 
is obviously extendable to many-particle Hamiltonians. It also enables one to generate new 
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QES Bose Hamiltonians by choosing another kinds of functions A, B, C, D in eqs. (p8|), 

The suggested approach can be useful in the problems of solid state physics when interac- 
tion between phonons is essential, quantum optics, theory of molecules, etc. In this context, 
especially important is the fact that our approach is extendable to many-particle systems. 
Concrete elaboration and applications of the obtained results deserve special treatment. 
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